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ABSTRACT 


The.ocgeanic response to time-dependent wind forcing approximately 100-200 km 
off the Northern California coast is examined through cross-spectrum analysis of 
current and wind time series. Current meter records 10 months in length from three 
moorings in deep water 100-200 km off the coast are analyzed in conjunction with 
surface wind analyses from the Fleet Numerical Oceanography Center (FNOC). 
Coherence is found between atmospheric and oceanic variables in the “locally forced” 
band (1-10 day period), in the “planetary wave” band (10-30 day period), and at low 
frequency (greater than 30 day period), though not in all of these bands for all records. 
The barotropic and first baroclinic dynamical modes appear to respond to wind forcing 
at different frequencies for two of the moorings analyzed. There is coherence between 
alongshore divergence and temperature fluctuations 100 km farther offshore, consistent 
with offshore advection by current filaments. Evidence of the Sverdrup balance is 
found for some periods greater than 18 days in the form of coherence between wind 
stress curl and the current component parallel to the local potential vorticity gradient. 
Complex bottom topography and the influence of coastal processes in the vicinity of 
the current meter moorings appear to greatly complicate the flow. There appears to be 
significant mesoscale variability in the region at scales too small to be resolved by the 
100 km spacing of the OPTOMA moorings. 
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I. INTRODUCTION 


A. BACKGROUND AND OBJECTIVES 

tie Navy has experienced a dramatic increase in the sophistication of its ship 
and aircraft based weapons systems and sensors over the last decade. Coincident with 
this increase in equipment capability has been an increased requirement to understand 
ocean structure and variability and their impact upon naval operations. The data 
currently used to describe and predict ocean structure are at best sparse, especially in 
light of the spatial and temporal scales typical of the upper ocean. An understanding of 
the relationship between atmospheric forcing and ocean currents is an important step 
toward describing and, ultimately, predicting such ocean variability with sparse oceanic 
and atmospheric data. 

The Ocean Prediction Through Observations, Modeling, and Analysis 
(OPTOMA) program is a joint Naval Postgraduate School/Harvard University project, 
sponsored by the Office of Naval Research, with goals of developing techniques for 
numerical forecasting of mesoscale ocean variability (such as eddies and meandering 
jets). OPTOMA has concentrated its efforts off the continental slope in an eastern 
boundary current regime. The present effort is focused upon a “test block” of ocean in 
the California Current System (CCS). OPTOMA is. developing an_ ocean 
descriptive-predictive system for four-dimensional data assimilation which could benefit 
directly from a better understanding of the relationship of ocean currents to 
atmospheric forcing. 

Previous studies have examined the oceanic response to atmospheric forcing in 
the mid-latitude North Pacific. Halliwell and Allen (1984) examined the response of sea 
level fluctuations in 1973 along the west coast of North America to alongshore wind 
stress at large scale (> 1000 km). They found atmospheric forcing to be most effective 
along northern California and Oregon. Forced fluctuations in sea level propagated 
poleward away from this forcing area, causing local sea level along the coast to be 
correlated with alongshore wind stress earlier in time and at a distant equatorward 
location. 

Niiler and Koblinsky (1985) examined long (~ 10 month) moored current meter 


records from far offshore in the eastern North Pacific with wind analyses from the U.S. 


Navy Fleet Numerical Oceanography Center (FNOC) and found wind stress curl to be 
significantly correlated with deep (below 150 m) ocean currents at periods of 10-100 
days. Much closer to the coast, Noble et.al. (1986) analyzed current meter records 
over the northern California slope and adjacent basin and found some weak evidence 
of wind-forced flow in both areas at periods less than 10 days. 

Stabeno and Smith (1986) analyzed current meter data from 11 moorings that 
were part of the Oregon State University/Sandia National Laboratory Low-Level 
Waste Ocean Disposal (LLWOD) Program (Heath eft. al., 1984) and 3 moorings that 
are part of the OPTOMA program. They found evidence of lens-like eddies in the 
currents within 1500 m of the bottom in the LLWOD records (100-200 km seaward of 
the 4000 m isobath), and much less low frequency (eddy) energy in the near-bottom 
OPTOMA records (shoreward of the 4000 m isobath). They also found evidence of the 
influence of CCS cold filaments originating near Point Arena in the near-surface 
OPTOMA records. 

The OPTOMA moorings analyzed by Stabeno and Smith (1986) were deployed in 
1984-85 off northern California by Oregon State University in association with the 
OPTOMA program (Smith et. al., 1986). FNOC surface wind analyses covering the 
mooring period and area are also available. The objective of this research is a detailed 
examination of the FNOC surface wind analyses and OPTOMA current meter data 
and the relationship between them. Specifically, the spectra of the current meter and 
wind records will be examined with the goal of better understanding how ocean 
currents respond over depth to wind forcing in the region equatorward of the 
Mendocino Escarpment. Knowledge of the oceanic response to atmospheric forcing 
from studies such as this will be useful in assessing and, perhaps, improving the 
modeling/prediction of wind-driven ocean circulation. The accurate prediction of ocean 
mesoscale variability (of which the wind-driven circulation is a component) 1s a vital 
prerequisite to the accurate forecasting of acoustic energy propagation that is needed 
by operational Navy antisubmarine warfare forces to efficiently and effectively carry 


out their missions. 


B. THE STUDY AREA 
1. Oceanography 
The current meter moorings are located within the OPTOMA NOCAL 


domain, an open ocean region nominally 300 km square. It 1s just offshore of the 


continental slope with a mean water depth of about 4000 m. The NOCAL domain is 
100-200 km south of the Mendocino Escarpment and is within 150 km of seamounts. 

The long-term climatological mean flow of the CCS is usually considered to be 
broad (~ 1000 km), shallow (~500 m), and slow (~4 cm/sec) (Hickey, 1978). The 
CCS forms the eastern boundary current of the North Pacific gyre and lies just to the 
east of the climatological North Pacific atmospheric ridge (see below), leading to the 
predomination of southward to southeastward winds during all but the winter season 
over the area south of 40N. These winds induce offshore surface Ekman transport 
and, consequently, coastal upwelling from February or March through August or 
September off Northern California. Due to a wind maximum ~ 300 km offshore they 
also produce a positive wind stress curl which induces Ekman suction (i.e., open ocean 
upwelling) from the coast to the offshore wind maximum. Similarly, there is a negative 
wind stress curl and open ocean downwelling seaward of the wind maximum. 

In contrast to its climatological mean flow, the instantaneous CCS is now 
known to be a very complex system of meandering jets, mesoscale eddies, and current 
filaments (Mooers and Robinson, 1984). The CCS coastal circulation structure varies 
seasonally, with northward flow in the winter, southward flow in the summer, and a 
northward undercurrent with a southward surface current in the summer and fall 
(Hickey, 1979). Meander patterns (which appear as large, cold filaments extending 
offshore) can be seen in satellite infrared images of the west coast of North America. 
From a series of dynamical model simulations, Ikeda and Emery (1984) concluded that 
these features stem from interaction of the mean current with bottom topography and 
grow because of baroclinic instability associated with vertical shear between the surface 
current and the undercurrents (during summer and fall). They further concluded that, 
in all seasons, irregularities in bottom topography generate meanders in the CCS with 
horizontal scales similar to that of the topography. There is also evidence from model 
results and hydrographic survey data for the existence of barotropic instability 
processes in the OPTOMA domain (Robinson ez. al., 1986). 

Philander and Yoon (1982) demonstrated, using numerical model experiments 
with both steady winds and fluctuating winds of various periods, that atmospheric 
forcing also plays a significant role in the complexity of eastern boundary current 
regions such as the CCS. Time dependent winds parallel to the coast create a vorticity 
source in the upper layers of the ocean near the coast which excite Rossby waves. 


These waves give rise to a complex system of poleward and equatorward currents near 


the coast that is superimposed on the mean flow (driven by the large-scale wind stress 
curl). 

Carton (1984) used the results of numerical experiments to show that oceanic 
response to isolated storms along an eastern boundary depends upon the ratio of storm 
duration to the time a coastal trapped wave takes to propagate to the point at which 
the response is measured. Carton found that the model ocean response to a storm was 
a vertically independent flow in the direction of the wind followed by an undercurrent 
flowing in the opposite direction. He suggested that variations in alongshore 
topography neglected in his study may be important sources of energy for 
vertically-varying waves. 

The combined effects of variability in the atmospheric forcing, the instabilities 
related to coastal effects (such as upwelling), and complex bottom topography make 
the nature of the ocean in this study area very complicated. 

2. Meteorology 

The properties of wind fluctuations along the west coast of North America 
vary over a wide range of time scales (Halliwell and Allen, 1986). Winter winds are 
primarily driven by propagating cyclones and anticyclones, and hence wind properties 
can fluctuate over time scales as short as several days. The summer winds are 
dominated by the effects of two quasi-permanent pressure systems - the North Pacific 
subtropical high and the southwest U.S. thermal low. Propagating surface cyclones 
(weaker than those of the winter) interact with the subtropical high to influence 
summertime winds. The subtropical high often strengthens and shifts toward the coast, 
causing intervals of strong equatorward (alongshore) winds off the California coast. 
The high occasionally builds far to the north for days to weeks at a time, preventing 


propagating cyclones from affecting the California coast. 
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Il. THE WIND FORCING 


A. WIND FIELDS 

Wind fields were extracted from the FNOC surface wind analyses for 16 grid 
points off the northern California coast, Figure 2.1, for the 13-month period June 1984 
through June 1985. For ease of reference, the grid points are designated A through P. 
The FNOC surface wind analysis and procedures used to construct the wind, wind 
stress, and wind stress curl time series used in this study are described in Appendix A. 

1. Seasonal Variability 

Monthly mean wind shows a distinct seasonal variability, Figure 2.2. In 

summer (represented by May 1985 in Figure 2.2), the mean has an equatorward 
component at all grid points. Within about 300 km of the coast the mean tends to be 
aligned with the coastline. The mean becomes more directly southward farther away 
from the coast. This anticyclonic mean flow pattern is associated with the North 
Pacific subtropical high (Halliwell and Allen, 1986). In virtually all cases, the mean is 
significant at the 95% level. The standard deviation ellipses are generally oriented with 
their major axes alongshore. The major axis direction, 9, is that of minimum 


Reynold’s stress, defined by 
yo Lt eae See ey 


where u’ and V’ are the perturbations from the time mean of the zonal and meridional 
wind components (respectively), and the angle brackets denote time-averaging. The 
standard deviations are generally larger for the more poleward grid points and are also 
larger close to the coast. The larger standard deviations observed at the poleward grid 
points reflect the position of the storm track which is poleward of the subtropical high: 
synoptic-scale atmospheric disturbances propagating through this region account for 
the larger observed standard deviations. 

In winter (represented by November 1984 in Figure 2.2), the mean 1s smaller 
and the standard deviations larger than in summer and frequently the means are not 
statistically significant at the 95% level. The winter storm track for synoptic-scale 
atmospheric disturbances is southward of its summer counterpart, with many major 


storms passing through the study area. Synoptic-scale fluctuations in the surface 


‘ t) 

« ° . ° e ty e to) * ° . 

. . . * e 8 s ° e . ° 

. r) ry r) ‘ ° . ° . . . 
42N ° . ‘ ‘ . . : . . . . 

i ae ee eee ee ee ee ee ee ee ee ee) 

® ® ° ° e ° t) ° ° a e 

* C) . . . * * . . . * 

) ty e ® ° a . . . s ° 

a s . ® e s A . . *. 


urek : 
ope Mendocino 


7] Ua sae toes cea aeaE EEE SEEEOSEREE WV LY MERESEEEEE SEED 
[ee et ee 


er ee ee 


ee ee ee ee ee ee eo 06 


SEN | 6 eee ecchtd eee 


sere Feces cede cece Feceee er Banc eeMoeecezabeseee Sees 0 +8 gmeoeece-ase8ereeeeessereterseee Je ge stu s + wows 


SAN ienteces a Ramco en ee 


ee ee ee eee Te eee ee ee Ur ee ee eee ee eee eee eee Cee ee 


see ee Peet eee es. cee Peees eee ses ee Pues eee Bee tae tee e es Bee = me =cewee eG eeese Peeeesrebesese Pu ee wearHeese= Hs sta scehuasewe Mee soe :Foeeee He oven 


slaieieis. 00.0 eu. o.g.m'e € 0: eleien wie'migiu'sio eects cle oecth cleteletierwieeifieln pies e/oieln/nlslensve/eielelsleteiel stele cieleieipisielsi=ielclniele'sinielgiaiaielsislstielatels(sipicla(eic(ole-s 6 sistoll'y 6 otatatelatet=lelert 





1SEW 1335WN 151W 129W I2Z7W ZEWieiZe Geil voi 


Figure 2.1. FNOC grid point and current meter mooring 
locations. [The dashed line outlines the domain 
of the OPTOMA]1] survey during summer 1984. 
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Figure 2.2 Monthly mean wind with standard deviation ellipses 
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for May 1985 (top) and November 1984 (bottom). 
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atmospheric pressure gradient caused by the propagation of storms through and near 
the study area alter the magnitude and direction of surface winds. These storms occur 
every few days in winter, accounting for the larger deviations and many changes in sign 
in winter wind velocities over periods of a few days, and, consequently, for the smaller 


mean observed in winter. 


B. WIND STRESS 

Wind stress is computed using the bulk aerodynamic method with the neutral 
flux drag coefficient of Large and Pond (1981), Appendix A. Wind stress is a more 
dynamically relevent variable than is wind and hence, all future references will be to 
wind stress unless otherwise noted. 

1. Seasonal Variability 

Largest magnitude events in the wind stress time series generally occur in 

winter and spring months, Figure 2.3. There are reversals in the sign of meridional 
component of the wind stress time series from October 1984 through March 1985. The 
meridional wind stress is equatorward from April through June 1985 with the exception 
of a few short, small magnitude poleward events, Figure 2.3. The zonal mean 1s 
eastward and the meridional mean is equatorward at all points, Table 1. The mean 
wind stress is different in direction near the coast from the general pattern over water, 
possibly indicating that the FNOC wind fields are sensitive to coastal topography, 
although Halliwell and Allen (1986) show that the winds derived from FNOC pressure 
fields do not correctly represent the coastal boundary layer effects. 

2. Autocorrelation Functions and Time Scales 

The wind stress time series are low-pass filtered using filter weights as 

suggested by Thompson (1983) with a cutoff frequency of about | cpd and half-power 
point at about 0.6 cpd (e.g., ~40 hour period). From complex autocorrelations 
(calculated for lags up to 68 days, 1.e., N/4, where N is the record length), the integral 
correlation time scale of wind stress for grid points G, H, K, and L are typically 1.5 - 
2.0 days for the zonal component and 2.5 - 3.0 days for the meridional component, 
Table 2. Wind stress at the point over land (point H) has an integral time scale of 
about 3.0 days for both components. Here, integral time scale, t; (as suggested by 


Richman et. al., 1976) is given by 


K 
a 2 
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TABLE 1 
WIND STRESS STATISTICS 


GRID-PT oe One 


0.3 0.6 
0.1 0.4 
0.3 0.8 
0.1 0.8 


* units: dyn/cm? 
< > denotes mean 


oO: standard deviation 


TABLE 2 
INTEGRAL TIME SCALES AND DECORRELATION TIMES 
FOR WIND STRESS FIELDS 
INT. TIME SCALE * DECORR. TIME * 


GRID PT ty 


2 
Pe 
Me 
oe 


*: units are days 





where R, is the correlation at lag kAt and K ts the maximum lag used in the 
correlation analysis (theoretically, a large lag at which the autocorrelation function has 
died out). The integral time scale is a measure of the time after which a new 
independent event may be considered to occur. The decorrelation time, Wap is the time 
for which 


R(Ty) = I/e. 
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These time scales indicate fluctuations in the wind records with 1-to-2 day time scales, 
consistent with the period of wintertime synoptic-scale atmospheric disturbances in the 
study area. 

3. Spectrum Analysis 

Characteristics of the wind stress field in the frequency domain are studied 
utilizing the rotary spectra techniques described by Gonella (1972) and Mooers (1973). 
The filtered wind stress time series at FNOC grid points G, H, K, and L are subdivided 
into five 360-point pieces with 50% overlap. Linear trends are removed from the 
series, a Hanning window smoother is applied to the spectrum of each piece, and piece 
averaging is performed to yield spectrum estimates with 10 degrees of freedom. 

The wind stress spectra at all four grid points are similar. All four spectra are 
white for periods greater than about 10 days and red for periods less than 5 days, 
Figure 2.4. The spectra of wind stress at points K and L have a peak at periods 
between 5 and 10 days, Figure 2.4. The energy level of the spectrum at point H (over 
land) is significantly lower at all frequencies than those of the spectra of the other three 
points, Figure 2.4. In all cases the clockwise (anticyclonic) spectrum is not 
significantly different from the anticlockwise (cyclonic) at the 90% significance level, 
Figure 2.5. Equal clockwise and anticlockwise components of the spectrum 1s an 
indication of rectilinear vice rotary motion. The spectra of all meridional wind stress 
components (as represented by point G) are about an order of magnitude more 


energetic than their zonal counterparts, Figure 2.6. 


Con WIND SiReSS:CURL 

Time series of wind stress curl are computed (see Appendix A) for a nearshore 
and an offshore case over a 12-month period to examine the spatial variability of such 
time series derived from FNOC wind data. The average wind stress curl (V x Tt) over 
the area denoted by grid points F, G, J, and K (with nominal center at 37N, 128W) 1s 
computed as the offshore case. The nearshore case is computed over the area denoted 
by points G, K, and L (with nominal center at 39N, 125W). The additional truncation 
error introduced by using a 3-point finite difference approximation to compute V x T at 
the nearshore position has negligible impact on the curl spectrum (see discussion, 
Appendix A). Monthly means of these curl time series were compared to those of 
Nelson (1977), who computed long term mean wind stress curl values from ship reports 


for 1° squares over the CCS. These results generally agree with Nelson’s (i.¢., same 
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Figure 2.4 Total spectra_of wind stress at FNOC grid 
points G, F ee 
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Figure 2.5 Auto-rotarv yn cetata of the wind stress field at 
grid point G (ty nical of all ocean ud points). Thick line: 
total spectrum, dashed line: anticlockwise, thin line: clockwise. 
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Figure 2.6 Spectra of meridional (a) and zonal (b) 
wind stress components, grid point G. 
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sign and order of magnitude in most cases), Table 3, usually positive along the coast 
and negative offshore. The agreement between the monthly means of this study and 
those of Nelson’s climatology is greater for the offshore case than for the nearshore 
case, Figure 2.7. The differences noted between the results of this study and Nelson’s 
climatology, however, do not necessarily require that one of them be in error. Possible 
explanations of the differences are: 


¢ Nelson’s climatology gives the long term mean (over a period > 100 years), and 
the resultant smoothing could decréase the mean values. 


e Nelson's climatology is based upon surface observations that are subject to many 
smaller scale atmospheric enomena_such as variation of the. marine 
atmospheric boundary layer depth. The FNOC winds used in this study reflect 
Se ic scale atmospheric disturbances and are generally. insensitive to the 

e 


smaller scale processes that affect in situ surface observations (Halliwell and 
Allen, 1986). 


e Calculation of wind stress curl on large grids such as_the FNOC grid can 
underestimate the values by as much as 50% (Saunders, 1976). 


e 1984-85 is not necessarily a typical year. 

In the auto-spectra of the two (offshore and nearshore) wind stress curl time 
series, there appears to be less energy in the offshore V x t as evidenced by the lower 
energy in its spectrum at low frequencies (periods <12 days), and at periods of about 5 
and 2 days, Figure 2.8. These observations, together with those for the wind stress 
field discussed above, point to the conclusion that the FNOC wind analyses are 
sensitive to coastal meteorological processes (associated with coastal topography, etc.) 
but do not necessarily represent them accurately (see Halliwell and Allen, 1986), and 
that wind stress curl computed from them contains spatial variability on scales of the 
FNOC grid spacing (~ 330 km at the latitude of this study). However, the V x t time 
series will not define scales of this magnitude; more likely, wavelengths the order 1000 
km will be resolved. A time series of wind stress curl over the OPTOMA domain is 
computed from wind stress at points G, K, and L for October 1984 to July 1985 (the 
interval of current meter data to be used in this study), Figure 4.3. The wind stress curl 
time series is filtered as described above for wind stress. The autospectrum of the wind 
stress curl has peaks centered at about 25, 12, 8, and 3 day periods, Figure 2.9. The 
spectrum falls off as ~ 3/4 (where @ is the frequency) for periods less than three 
days, which is consistent with the spectral character of V x t selected by Frankignoul 
and Muller (1979) for model studies after examination of available wind and wind 


stress curl time series. 
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Figure 2.7 Comparison of wind stress curl monthly means 
for nearshore (top) and offshore (bottom) cases. 
SoOlicdime; presemms tar maashed IaeaNelson (1977). 
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re 2.9 Autospectrum of wind stress curl over 
the OPTOMA domain for the period covering 
the current meter moorings. 


Ii. THE CURRENTS 


A. CURRENT DATA 

The three OPTOMA current meter moorings were deployed off northern 
California from September 1984 to July 1985 by Prof: R.L. Smith at Oregon State 
University in association with the OPTOMA program. The moorings (designated M1, 
M2, and M3), Figure 3.1, were separated by about 100 km. Five Aanderaa current 
meters were installed on each mooring to record temperature, pressure, current speed 
and direction, at 1 hourly intervals. The deepest instrument at each mooring was 
approximately 200 m above the bottom. Calibration and post-processing were carried 
out at Oregon State University (Smith et. al., 1986). Two meters each on moorings M1 
and M3 failed or provided incomplete records. Hence, a total of 11 records from the 
three moorings are used for this analysis, Table 4. 

The most comprehensive time overlap of the current records is 0400, 02 October 
1984 to 2300, 02 July 1985 (inclusive) and, hence, this period is selected for the focus of 
study. All current records are low-pass filtered and then subsampled at six-hourly 
intervals. None of the current time series appear to be rotary in nature, Figures 3.2, 
3.3, and 3.4. Means of all records are not statistically significantly different from zero 


at the 95% level, Table 5. The 95% significance level is given by the standard error: 
26/N1/ 2 


where © is the standard deviation and N is the number of degrees of freedom. 
Standard deviation ellipses, aligned so that the major axis is parallel with the direction 
of minimum Reynold’s stress, show the maximum variance to be polarized in the 
north/south direction. Monthly standard deviation exceeds monthly mean in all cases. 
There are large month-to-month variations in the mean direction as well as magnitude, 
but in all except a few cases the monthly mean 1s not significantly different from zero 
at the 95% level, Figure 3.5. 
1. Autocorrelation Functions and Time Scales 
A prominent signal with a period of ~70 days dominates the autocorrelation 


functions of the meridional components of M2-340 and M3-350, Figure 3.6, as for 
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TABLE 4 


CURRENT METER DEPTHS AND RECORDED VARIABLES 
USED IN ANALYSIS 


DEPTH (m) VARIABLES 
MOORING M1 
MOORING M2 


MOORING M3 


LEGEND 


: Zonal velocity component 

: Meridional velocity component 

: Temperature 

: Denotes instrument ~ 200 m above bottom 





many other levels except those near the bottom.! The autocorrelation function of the 
zonal component of M3-350 has a markedly longer dominant period than do those for 
M 1-375 and M2-340, Figure 3.6. Decorrelation times for both components range from 
2-to-22 days; smallest values are for the near-bottom records at M1 and M2 (Table 6). 
Integral time scales range from 4-to-38 days. Again, the smallest values are for the 
near-bottom records at Ml and M2, and for the u-components of MI1-175 and 
M 2-145. 


1M1-375 denotes the current record from M1 at 375 m depth, Otherenni 
meter records will be referred to in a similar manner. 
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Figure 3.4 Current time series from mooring M3. Depths from 
top are 350, 800, and 1185 m. Units: cm/s. 
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TABEES 


STATISTICS OF CURRENT AND TEMPERATURE RECORDS 





RECORD =U 35 On <y>” Gc. 8 # <T> ** o7** 











M1-175 -2.4 4.8 -2.5 6.7 8 8.5 0.47 
M1-375 -0.9 ag -1.2 2.9 >> 6.3 0.26 
M 1-3250 0.2 les 0.4 1.5 -27 1.5 0.06 
M2-145 -1.4 Dal 1.0 8.9 -6 8.4 0.47 
M 2-340 -1.0 ZS has 5.4 -9 6.2 0.35 
M 2-800 -0.5 1.8 1.1 2.6 -22 4.2 0.11 
M2-1190 0.0 1.4 0.6 ea -34 a 0.05 
M 2-3560 0.0 0.9 -0.3 1.1 10 Fe: 0.01 
M 3-350 -3.5 4.6 -0.6 4.4 -59 6.3 0.38 
M 3-800 -2.4 3.4 -0.1 Pe -74 4.2 0.13 
M3-1185 -1.7 2.6 -0.8 2,5 -52 Sy 0.08 


< >: denotes mean 

o: denotes standard deviation 
* >) Um tS: ems 

% Seal TVGG ee 


#: in degrees clockwise from north 


2. Spectrum Analysis 
The current records are processed and transformed for spectrum analysis as 
described for the wind records in Chapter II. The auto-rotary spectra of all the current 
records except the near-bottom records are similar. The spectra do not display 
distinctive rotary character, except at the lowest frequencies and some isolated higher 


frequencies (®~0.22 cpd). At the lowest frequencies the clockwise component 
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Figure 3.5 Monthly mean and standard error of selected 
current records at \il, \l2 and N{3. Units: cmis. 
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Figure 3.6 Autocorrelation functions of zonal and meridional 
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TABLE 6 


INTEGRAL TIMES SCALES AND DECOR 
TIMES FROM CURRENT RECOR 


RELATION 
DS 


_ INT. TIME SCALE * DECORR. TIME * 
DEPTH (m) 
MOORING M1 
175 
375 
3250 


MOORING M2 


*: units are days 





dominates; at ®@ ~0.22 cpd the anticlockwise component is more energetic. The spectra 
are red to ® ~0.3 cpd and essentially white at higher frequencies, Figure 3.7. 

The spectrum of M1-3250 is white for @ <0.05 cpd and m>0.15 cpd. There 1s 
a large amplitude peak in the spectrum (both rotary components) centered at w~0.1 
cpd, Figure 3.8. The M2-3560 spectrum is essentially white to about 0.1, decreases 
sharply by about an order of magnitude, and is white at W>0.2 cpd, Figure 3.9. 

3. Dynamical Mode Decomposition 

The current field dynamical modes are shown for completeness, Figure 3.10, 
since they are used in analysis of the current response to wind forcing in Chapter IV. 
The current field is decomposed into dynamical modes using a Brunt-Vaisala frequency 
profile, N*(z), computed from an average of OPTOMA CTD casts within ~ 50 km of 


the moorings (Rienecker et. al., 1986). 


39 


= 


10 


(cm s_')’/cpd 





ioe 10 10° 
frequency (cpd) 


Figure 3.7. Auto-rotarv spectra of current at \f2-340, 
(tvpical of other records). Thick line: total 
spectrum, dashed line: anticlockwise, thin line: clockwise. 
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Figure 3.8 Asin Figure 3.7, but for M1-3250. 
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As in Figure 3.7, but for M2-3560. 
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B. SPATIAL VARIABILITY OF CURRENTS 
l. Vertical Coherence 

Planar cross-spectra of all the pairs of meridional current components at 
mooring M2 are computed to investigate the vertical coherence of the flow. The 
coherence-of each pair at a period of 45 days is computed and compared to the 
coherence of the other pairs as a function of vertical separation (Az), Figure 3.11. At 
low frequency, the coherence generally decreases with increasing Az at mooring M2, 
Figure 3.11. This variability, also found at MI and M3, contrasts sharply with the 
nearly barotropic flow found by Niiler and Koblinsky (1985) much farther offshore and 
north of the Mendocino Escarpment (42N, 152W), and suggests a much more 
complicated flow structure in the OPTOMA domain than for “open ocean” regimes. 

2. Horizontal Coherence 

The only meters at a common depth between the three moorings that are 
available to investigate horizontal variability are those at about 350 m. The planar 
cross-spectra of the principal axis components are computed as a measure of the 
horizontal coherence of the flow. MI and M2 are coherent” only at 8-9 day periods, 
and M1 is coherent with M3 only at very short (<3 day) periods. M2 and M3 are 
coherent for periods longer than ~40 days, Figure 3.12. Hence, there is limited 
coherence between all three moorings at the 350 m level, but at different periods for 
each pair. The moorings are located in an area of complex bottom topography, Figure 
3.1, which presumably accounts for the general lack of horizontal coherence in the 


current field. 


“Unless otherwise indicated, the 90% significance level will be used throughout 
for coherence/no coherence determinations. 
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Figure 3.11 Coherence between pairs of meridional current 
components at \12 as a function of vertical separation (Az). 
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Figure 3.12 Coherence between current components parallel 
to the principal axis direction of M1-375, M2-340, and 413-350. 
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IV. THE CURRENT RESPONSE TO WINDS 


Theoretical studies and model results (e.g., Muller and Frankignoul, 1981 and 
Willebrand et. al., 1980) suggest that the relationship between atmospheric forcing and 
the open oceanic response is complex and not readily observable. Muller and 
Frankignoul (1981) compared model results to observations for an open ocean (weak 
eddy field) regime. Their results suggest that atmospheric forcing may be the dominant 
energy source for geostrophic eddies in regions removed from strong currents. They 
concluded, however, that direct evidence of this may be difficult to establish, especially 
at low frequency. Willebrand ez. al. (1980) used wind data derived from synoptic 
NMC surface pressure charts to force a closed-basin numerical model. From the 
barotropic vorticity equation for a flat-bottom ocean, they concluded that the oceanic 
response to stochastic atmospheric forcing 1s characterized by three distinct frequency 
ranges: 

e At periods of 1 to 10 days the response 1s locally forced. 

e At periods of 10 to 30 days the response 1s 1n the form of planetary waves. 

e At periods longer than a month the motion corresponds to a Sverdrup balance. 
Their results from a nonlinear numerical model forced by realistic winds confirmed this 
conclusion. They also found that bottom topography significantly altered coherence 
between the ocean and atmospheric forcing. Their conclusion that the Sverdrup 
balance dominates only at periods greater than 30 days contrasts with the results of 
Niiler and Koblinsky (1985), who found evidence of a time-dependent Sverdrup balance 
at periods of 10-100 days. 

The cross-spectra of wind and current meter records described in Chapters I] and 
II] are analyzed seeking: 


¢ coherence between ocean currents and wind stress in each of the frequency bands 
defined by Willebrand ez. al. (1980 


e evidence of Ekman pumping 
e evidence of the Sverdrup balance 


e other factors that may selectively mask the oceanic response to wind forcing in 
frequency, space, or time. 


It is anticipated that this regime is different from those previously investigated. It 1s 


near (100-200 km) the coast, yet probably too far offshore to observe coastal eflects 
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such as coastally trapped waves. On the other hand, it is clearly not analogous to the 
Open ocean regimes studied by numerous earlier investigators. The bottom topography 
in the study area is very complicated, Figure 3.1. It is anticipated that this will tend to 


destroy (or cover) coherence between the ocean and wind forcing. 


A. WIND STRESS/CURRENT CROSS-SPECTRUM ANALYSIS 

For consistency, wind stress at point G will be used for all wind/current 
cross-spectrum analyses using the OPTOMA current data. Point G is closest to all of 
the mooring locations, Figure 2.1. Standard deviation ellipses of points K and L are 
larger than that of point G, Figure 4.1, reflecting their position closer to the storm 
track across Oregon and northern California. Since the moorings are somewhat south 
of the storm track, point G wind stress is probably more representative of the actual 
wind over the moorings. The orientation of the standard deviation ellipse for point G is 
steadier with frequency than those for points K and L and is aligned ~15° west of 
north, Figure 4.1. This alignment is consistent with the wind direction that is expected 
on the downstream side of the North Pacific subtropical high (see Chapter I). The 
standard deviation ellipses of points K and L are aligned east of north and, hence, not 
consistent with the wind flow expected around the subtropical high, Figure 4.1. 
Cross-spectra of wind stress at points G, K, and L with several current records show 
point G to be the most coherent with current records. In several ways, then, wind 
stress at point G appears to be the most representative of wind stress over the 
OPTOMA moorings and 1s selected over points K and L. All further references to wind 
stress will be to wind stress at point G unless otherwise noted. 

The wind stress principal axis direction is aligned about 15° west of north and 1s 
nearly constant with frequency, especially over portions of the spectrum where ellipse 
stability is high. Neither the wind stress nor any of the current time series 1s 
significantly rotary in nature, so the planar cross-spectra of Tt, vs. u and ty vs. Vv are 
computed. The planar cross-spectra have many coherent peaks, but more for the ty 
vs. v cases than for the t, cases (probably since ty is more energetic than t,), Table 7. 
There are coherent peaks in the “locally forced” band (periods of 1-10 days) in each of 


these cross-spectra, except T,, vs. v at M2-800. These peaks are concentrated at periods 


y/ 
of 3-9 days, Table 7. There is very little coherence in the “planetary wave” band (10-30 
day periods); M1-175 and M1-3250 for the Ty vs. V case and M2-145 for the t, vs. u 


case are the only ones with coherence in this band. There is coherence between wind 
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Figure 4.1 Mean vectors and standard deviation ellipses for 


wind stress at FNOC grid points A through P over the 
period of the study. 
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H7W 


stress and current components in the “Sverdrup balance” band (<30 day periods) for 


several of the cases, Table 7. 


TABLE 7 
COHERENT 


PEAKS BETWEEN WIND STRESS AND 
CURRENT COMPONENTS 


MOORING/DEPTH (m) PERIODS OF COHERENT PEAKS (days) 


Ty Vs. U ty VS. V 


M 1-175 | 7* 5.5, 10, 15, 22 
M 1-375 3.5*, 7% 79 
M 1-3250 4,7 5*, 6.5*, 9, 15 


M2-145. 4*, 8*, 24* 3.3, 5*, 9, > 45 
M2-340 15 5,9, >45 
M2-800 3.3 43-7] 
M2-1190 5 9, 43-55 
M2-3560 6* 


M 3-350 a Sees 
M 3-800 3. SEA 345 
M3-1185 Lee) 


*: coherence exceeds 95% level 





The current components parallel to the principal axis direction for M1-375, 
M 2-340, and M3-350 are designated M1-375 4 33, M2-340 19, and M3-350_¢69, where 
the subscript denotes the principal axis direction in degrees clockwise from north. 
M1-375 4.33 is coherent with t, ata period of ~8 days. M2-340_,9 has coherent 


peaks with T,, at 5.5 and 9 day periods, and for all periods greater than 45 days. There 


M 


is significant coherence between t, and M3-350_¢9, especially at a period of 9 days, 


y 
Table 8. In general, computation of cross-spectra using current and wind stress 


components parallel with their respective principal axis directions does not appreciably 
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affect the character of observed coherences from those in Table 7. As with the 
zonal/meridional cross-spectrum analysis described above, there is evidence of ocean 
response to wind forcing at periods of 1-10 days, at periods of 10-30 days, and at long 
(> 30 days) periods. Higher frequency coherence (periods <10 days) is consistently 
present in cross-spectra computed, intermediate frequency coherence is present in only 
one case, and lower frequency coherence (period > 30 days) is present in only one case. 
This is consistent with the results of Muller and Frankignoul (1981), who expected 
ocean/wind coherence to be small and difficult to observe in the more energetic lower 
frequencies because of destructive interference of different wavenumbers, and with 
Willebrand et. al. (1980), who expected to find no coherence in the 10-30 day period 
band. 


TABLE 8 
COHERENT PEAK 
PARALLE 


S BETWEEN WIND STRES 
L TO PRINCIPAL AXIS DI 


S 
RECTIO 


MOORING/DEPTH (m) PERIODS OF COHERENT PEAKS (days) 
M1-375 4 33 4.3, 7.5*, 10, 15 
M2-340_19 5.5,9, > 45 
M3-350_¢9 65 9.5.5, 10% 


*. coherence exceeds 95% level 





B. CROSS-SPECTRUM ANALYSIS WITH WIND STRESS CURL 
1. Wind Stress Curl Versus Temperature and its Time Derivative 

The cross-spectrum of wind stress curl with temperature records 1s computed, 
seeking evidence of Ekman pumping. V xT is coherent with temperature at M1-375 
for 14 and >67 day periods, and with temperature at M2-340 for 18 and >6/7 day 
periods. Temperature at M3-350 is coherent with V x t for 6 and 28 day periods (the 
28 day peak is significant at the 95% level). Essentially no coherence is found between 
V x t and temperature records at any other level, but the coherence found at the 350 m 


level may be evidence of the Ekman pumping mechanism. 


The heat equation may be expressed as: 


a = pes -ueVT + K*v?T 1 ae 
t OZ pc, 
(1) (2) (3) (4) (5) 


where term (1) is the local time derivative of temperature, (2) is vertical advection , (3) 
is horizontal advection, (4) is diffusion, and (5) is diabatic heating. If terms (3) through 
(5) can be assumed small and ignored, the local time derivative of temperature (dT/ dt) 
may be taken as a measure of vertical velocity, w, which is related to wind stress curl 
by 


w= Vx (t/O, 


where p is the water density, f is the coriolis parameter, and Tt is the wind stress. 
Usually the spatial variation of wind stress is greater than the variation of the coriolis 


parameter with latitude, so vertical velocity may be approximated as 
w = (1/pf) V x f. 


Thus, the Ekman pumping vertical velocity is proportional to the wind stress curl (Gill, 
1982). 

The cross-spectrum of wind stress curl and 0T/ét derived from the ~350 m 
temperature records at all three moorings, and from the shallow and deep temperature 
records at M1 and M2 are computed. V x Tt is coherent with @T/ét at M3-350 for 
periods of ~30 days and 8-9 days. At M1-375 and M2-340 there is no significant 
coherence with V x t for any period. V x t has a coherent peak with dT/dt at M1-3250 
for 14-day period. There are also a few peaks at high frequency (<3 day period) for 
M1-175, M1-3250, M2-145, and M2-3560. The general lack of coherence between 
wind stress curl and 0T/dt suggests the importance of the neglected terms in the heat 
equation, especially (in this region) horizontal advection. 

2. Wind Stress Curl and Temperature Versus Current Horizontal Divergence 

Horizontal divergence is computed from the 350 m current at all three 
moorings. The current horizontal divergence spectrum is red at low frequency, has a 


peak at 10-day period, and becomes white for periods <5 days. 
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Wind stress curl and horizontal divergence are coherent at periods of about 5 
and 8 days. There is another smaller peak at high frequency (in the “weather” band), 
Figure 4.2, but in general there is little coherence between V xt and horizontal 
divergence. The V x T time Series is much noisier than the horizontal divergence time 
series, and no correlation between the two is readily apparent, Figure 4.3. Similarly, 
there is very little coherence between temperature at any of the ~350 m records and 
horizontal divergence, and no apparent correlation between any of them and horizontal 
divergence, Figure 4.3. 

Since the expected coherence between horizontal divergence and temperature 
is not found, an alternate hypothesis is proposed. If alongshore convergence between 
M1 and M2 at 350 m were to be accompanied by offshore flow (toward M3) of the 
same magnitude, then the horizontal divergence computed between M1, M2, and M3 
would be small. Thus, significant alongshore convergence/divergence would necessarily 
not be reflected in the total horizontal divergence. There are several alongshore 
convergence events that are followed in time by a cooling event at M3-350 (~ 100 km 
farther offshore), Figure 4.3, and several examples of the opposite (alongshore 
divergence leading a warming event offshore). Alongshore divergence and temperature 
at M3-350 are coherent for periods greater than about 26 days, Figure 4.4. Their phase 
(with alongshore divergence leading temperature at M3-350) increases with increasing 
period (at periods for which they are coherent) and is greater than 180°, which is 
consistent with advection in the offshore direction. This is an indication that advection 
is an important process affecting variability in the study area, at least in the upper 350 
m of the water column. This coherence is consistent with the CCS current filaments 
extending offshore, and to depths of 350 m that were described by Mooers and 
Robinson (1984) and Ikeda and Emery (1984). 

3. The Sverdrup Balance 
The Sverdrup balance in the presence of bottom topography 1s given by 


Ue WAH) = (pHy! V xt, 


where U is the horizontal current vector, H is the depth of the ocean, and V 1s the 
two-dimensional horizontal gradient operator. This balance relates the current 
component across potential vorticity (f/H) contours to local wind stress curl. To seek 
evidence of this balance, the cross-spectrum of wind stress curl with 350 m current 


components varied by 30° are computed to determine the direction of maximum 
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Figure 4.2. Coherence and phase between wind stress curl 
and horizontal current divergence at 350 m. 
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coherence for each mooring. For M1-375, 30° east of north is the direction of 
maximum coherence - not inconsistent with the direction of the local bathymetric 
gradient, Figure 3.1. Wind stress curl is coherent with current parallel to V(f/H) at M1 
for periods of 18 and 28-50 days, Figure 4.5. Wind stress curl and M2-340 (with 
current component 120° east of north) are coherent only at 90 days - this is not strong 
evidence for a Sverdrup balance, Figure 4.6. Wind stress curl is incoherent with 


M3-350 current components of any orientation. 


C. WIND STRESS/DYNAMICAL MODE CROSS-SPECTRUM ANALYSIS 
Component-wise cross-spectra of wind stress with current dynamical mode 
amplitudes are computed to investigate separately the barotropic and baroclinic ocean 
response to wind forcing. 
1. Mooring M1! 
The meridional component of the barotropic mode is coherent with ty at 
periods of <3, 5, 6.5 and 10 days, Figure 4.7. The v-component of the first baroclinic 


mode has a large coherent peak with T,, at a period of 15 days and a small peaks at 6.5 


y 
and 3-4 days, Figure 4.8. Hence, both modes share a coherent peak with t,, at 3-4 


y 
days and at 6.5-day period. At lower frequencies, both modes are coherent with the 
wind forcing, but at decidedly different frequencies. 
2. Mooring M2 


The v-component of the barotropic mode and t,, have coherent peaks at 


isolated periods <5 days and centered at 8 days. They are Sis coherent over periods 
greater than 55 days (just barely at the 90% level), Figure 4.9. The v-component of the 
first baroclinic mode is coherent with ty at periods <5 days and > 55 days. It also has 
a coherent peak centered at 5 days, Figure 4.10. Hence, both dynamical modes are 
coherent with ty at short and long periods, but respond to f,, at different intermediate 


uy 
frequencies. These intermediate frequencies are different from those at M1. 


i DISCUSSION 

Zonal and meridional wind stress components are coherent with their 
corresponding current components at various periods in all three frequency bands 
defined by Willebrand er. al. (1980). The same is true of wind and current components 
(at 350 m depth) aligned with their respective principal axis directions. 

There is weak evidence of Ekman pumping in the form of coherence between 


wind stress curl and temperature at the ~350 m level, and essentially the same 
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Figure 4.5 Coherence and phase between wind stress curl 
and current component 30° east of north at M1-375. 
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Figure 4.6 As in Figure 4.5, but for i2-340 and 
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Figure 4.7. Coherence and phase between meridional wind stress 
and meridional component of barotropic mode, M1. 
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Figure 4.9 Asin Figure 4.7, but with the 
barotropic mode, M2. 
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between OT/Ot and V x t. Neither wind stress curl nor temperature is coherent with 
horizontal divergence computed at the 350 m level, but alongshore divergence between 
M1 and M2 is coherent with temperature at M3-350 (farther offshore) at low 
frequency. Horizontal advection masks the evidence for Ekman pumping. To see if 
masking of the current response to wind stress curl may be due to coastal effects, the 
results are compared to those for mooring CMMW10, one of the LLWOD moorings 
analyzed by Stabeno and Smith (1986), Figure 3.1. CMMW410, with five current 
meters, was deployed approximately 200 km northwest of the OPTOMA mooring 
locations for ~10 months in 1982-83. There is very little coherence between V x t and 
either temperature or OT/Ot at mooring CMMW1O0, however, CMMW1O0 was located 
very close to the V xt longitude off the northern California coast (Nelson, 1977). 
Hence, the much less energetic V x t forcing over CMMW10 may account for the lack 
of coherence with temperature and @T/ét. 

At moorings M1 and M2, the barotropic and first baroclinic dynamical mode 
components are coherent with corresponding wind stress components, but at different 
periods for the two modes. The v-component of the barotropic mode at CMMWI1O0 is 
coherent with ty at periods of 8 and 16 days, Figure 4.11. The v-component of the first 
baroclinic mode at CMMW10 is coherent with ty at the same periods, Figure 4.12, in 
contrast to M1 and M2 where the barotropic and baroclinic responses were at different 
periods. This could be a result of simpler dynamics in the CMMW10 regime (simpler 
bottom topography, farther from the coast) than in the OPTOMA domain. For 
example, the ratio of long-term integrated energy in the barotropic and first baroclinic 
modes is 0.64 for M2 and 1.84 for CMMW10. However, M2 and CMMW1O0 were 


deployed in different years, and this alone could account for the observed differences. 
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Figure 4.11 Asin Figure 4.7, but with the 
barotropic mode, CMMWI1O0. 
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V. SUMMARY AND CONCLUSIONS 


A. SUMMARY 

This study has examined the oceanic response to wind forcing by analyzing 
moored current meter records with FNOC surface wind analyses. There is coherence, 
to one degree or another, between wind forcing and ocean current variables in all three 
frequency ranges defined by Willebrand et. al., (1980), and the results are summarized 
by these ranges. 

1. Locally Forced Response (1-10 day period) 

There are many coherent peaks in the cross-spectra of wind stress and ocean 
currents. These peaks are concentrated at periods of 3-to-9 days and are very robust, 
being in the component-wise planar cross-spectra of wind stress and current for every 
mooring and level. The large grid spacing of the FNOC wind analyses (~330 km at 
the latitude of this study) is probably too large to resolve all of the wind features that 
are forcing the oceanic response in this frequency band, which may account for some 
of the imperfection in wind stress/current coherence in this band. The importance of 
horizontal advective processes in the study area may also account for some of the 
imperfections in wind stress/current coherence in this band, although, 1n this study, the 
evidence of the advective processes is at much longer periods. 

2. Planetary Wave Response (10-30 day periods) 

Very little coherence between atmospheric and oceanic variables is found in 
this range, except for M1-3250. The autospectrum of M1-3250 has a large peak 
centered at a period of ~13 days, which may represent the free, bottom-trapped 
topographic waves observed in the model results of Willebrand er. al. (1980). These 
bottom-trapped topographic waves were generated in their model by vertical motion 
produced as the wind-induced flow interacted with topography. The large peak in the 
M 1-3250 spectrum in this frequency range is possibly evidence of this process. 

Limited evidence of Ekman pumping ts also found in this frequency range. 
Temperature changes at 350 m depth are coherent with wind stress at periods of 14, 18, 


and 28 days, as well as at long periods (> 67 days). 
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3. Sverdrup Balance Response (> 30 day periods) 

Evidence of the Sverdrup balance is found at MI, where U ¢ V(f/H) is 
coherent with wind stress curl for periods of 18 and 28-50 days. The M1 results agree 
with those of Niiler and Koblinsky (1985), who found evidence of the Sverdrup balance 
at 10-100-day periods (albeit much farther offshore). The coherence between wind 
stress curl and U e V(f/H) for other periods at M1 and for all periods at M2 and M3 is 
most likely destroyed by the horizontal advective processes that appear to dominate 
variability in the upper levels of the water column in the study area. 

4. cffects of Bottom Topography 

The degree of coherence between wind forcing and current variables varies 
markedly horizontally and vertically. Model results suggest that bottom topography 
may account for much of this inhomogeneity. The results of numerical experiments 
with smoothed bottom topography, homogeneous ocean, and stochastic wind forcing 
show considerable horizontal variation in the current spectrum over short distances 
(Willebrand ez. al., 1980). Similar experiments with a stratified ocean demonstrated that 
topography leads to mode coupling and energy exchange between modes (Muller and 
Frankignoul, 1981). It is likely that the fairly complicated bottom topography in the 
vicinity of the OPTOMA moorings significantly alters the flow and partially destroys 
any coherence with wind forcing. 

5. Coastal Effects 

The proximity of the coast to the OPTOMA mooring locations is a major 
complicating factor that may partially destroy wind-current coherence. Satellite infrared 
images of the Oregon and northern California coast contain many cool filaments and 
meandering jets extending offshore through the OPTOMA domain (Ikeda and Emery, 
1984). These cool filaments and meandering jets have accompanying current filaments 
below the surface which can also be seen in hydrographic data (Mooers and Robinson, 
1984). Alongshore divergence between M1 and M2 at 350 m depth is coherent with 
temperature changes at M3 (100 km farther offshore), consistent with offshore 
advection by CCS current filaments. These results indicate that the CCS current 
filaments extend at least as deep as 350 m. 

In summary, the oceanic response to wind forcing in the OPTOMA domain ts 
very complex. The results of this study are not inconsistent with earlier model and 
observational studies. The close proximity of the continental slope (about one Rossby 


radius away) and the strong advective processes that appear to dominate upper level 
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variability in the study area make it definitely unlike a mid-ocean regime, but it is too 
far removed from the coast to be affected by coastal jets as observed by Kundu and 
Allen (1976), coastally trapped waves, or other strictly coastal processes. Hence, the 
OPTOMA domain is neither mid-ocean nor coastal in nature. 

Overall, the issues raised here are important for the development of ocean 
circulation models, observing systems, and ocean prediction systems for eastern 
boundary current regimes. To date, the OPTOMA program has been able to ignore 
atmospheric forcing and to proceed without current meter data in ocean prediction 
studies. This has been at the expense of using much XBT and CTD data. The present 
results provide the hope that an atmospheric forcing prescription could possibly be 
economically traded-off, with the use of appropriate models, against mass field data in 
future ocean prediction studies. Presumably, there would also be an important role, yet 


to be established, for current meter data in such studies. 


B. RECOMMENDATIONS 

This study is a first step in describing the oceanic response to wind forcing in an 
eastern boundary current regime like the OPTOMA domain. Future studies should 
examine the OPTOMA current data for coherence with distant winds, as found in 
model studies by Carton (1984) and in observational studies of wind stress and sea 
level response very close to the coast by Halliwell and Allen (1984). Analysis of the 
OPTOMA current records should also be undertaken with wind forcing information 
from other available sources (such as meteorological buoy wind data and the National 
Meteorological Center’s LFM and NGM wind analyses) that may resolve wind 
features of smaller scale than the scales resolvable by the FNOC grid spacing. The 
results of this study show that there is a considerable amount of mesoscale ocean 
variability in the OPTOMA domain that is of too small a scale to be resolvable by the 
100 km spacing used for moorings M1, M2, and M3. Clearly, more (and more closely 
spaced) current meter moorings will be required to describe fully all of the horizontal 
variability in this region. Current records closer to the surface will be required to 


adequately sample the effects from the coastal transition zone. 
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APPENDIX A 
PROCESSING AND ANALYSIS OF FNOC WIND DATA 


FNOC surface wind analyses, on a 63x63 Northern Hemisphere polar 
stereographic grid, are generated at 6 hour intervals using NOGAPS (Navy 
Operational Global Atmospheric Prediction System) forecast fields as first-guess 
analysis fields, to which new data are added using a form of the successive corrections 
interpolation technique (Rosmond, 1981). These analyses are archived on a 6 hour 
cycle (00Z, 06Z, 12Z, 18Z) and are used to re-initialize the NOGAPS global 
atmospheric forecast model. 

The time gaps in the FNOC analyses were identified and the duplicate records 
deleted on a first editing pass. Data points that were clearly spurious in magnitude 
(differed from the mean by more than 30, where o is the standard deviation) were also 
omitted. Data gaps were found to be short (of 24 hours or less), and they were filled 
by simple linear interpolation. The use of a higher order polynomial or statistical 
interpolation would probably be required to fill larger gaps satisfactorily. 

An iterative technique which assumed a neutral flux profile was used to adjust 
the FNOC values, which are at a height of 19.5 m, to 10 m, at which the wind stress 
was computed. No change was made to wind direction. The computation used the 
neutral flux drag coefficient of Large and Pond (1981), and the iterative scheme is 


outlined below. Given: 
Uj;g: u wind component from FNOC data 
V9: V wind component from FNOC data 
K: Von Karman’s constant 
€: convergence criteria for adjustment factor 
Cy: drag coefficient 
ux: friction velocity 
g: acceleration due to gravity 


Let: 
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K = 0.4 
&€= 0.1 
g = 9.8 m/s 


Begin the iteration: 


(10Cq)y, = 1.14, |Vlp < 10 m/s 


0.49+0.065|VI,, |Vi, 210 m/s 


(Us) = (us) (Cg)l/? 


[In(10.0g/ux)/K] - 10.6 


1 “Wn(9.5g/ux2)/K] - 10.6 
Vig = [Von 


The iteration continues until VVIn-TVi,_ pl < € Then, the adjustment factor, pt, is 


given by: 
ee IVEY | 6 
Finally, 
“one U9 
Ores & * 19> 


The FNOC polar stereographic projection 1s true at 60N. Using indicies 0-62, 
the pole is at (31,31). The 80W meridian 1s aligned with 1=31 on the grid. The grid 
spacing, Ax, 1s 381 km at 60N and varies with latitude: 


Ax = (381)[(1+sin 9)/1.866], 


where @ 1s the latitude. Given the (1,j) coordinates of a polar stereographic grid point, 
latitude and longitude (A) of the point may then be computed by (e.g., Haltiner and 
Williams, 1980): 

d 


10° + tan” (yilx) 


@ = sin ![(z*-y,*)/(2"+y,)], where 


7] 


a{1+(sin 60°)]sin 0, 


N 
| 


0m) 
I 


Oe 

a = 6371 km (radius of a sphere with same volume as earth), 
xX = (i-31)381, and 

i (j-31)381. 


The u and v wind components (Un gVps) obtained from FNOC data are referenced to 
the polar stereographic grid coordinate frame. The u and v components in geographic 
(1.e., E-W/N-S) coordinates (u_,, V,,) are given by: 

v= “Up,(sin y) + Vp s(CO8 Y) 

ye “Up,(Cos Y) - Vps(sin Y) 
where Y = tan”!(y,/x}). Wind speed and direction information that are relative to true 
north are also available from FNOC. 

The processing outlined in this appendix was applied to the FNOC surface 
analysis data to generate a wind time Series at 10 m for each grid point of interest. The 
wind stress may now be computed: 

t, = p,CqUujo, and 

ty = pP,C qgUV19, where 
\7 m2 es 10°9 gm. cm”? (air density), and 
WW (U49° ar vig, 


and uyjqg and V)¢ are in either geographic or polar stereographic coordinates. 


Haney ez. al. (1981) gives the formula for the wind stress curl: 
Vxt= m*[4/8x(ty/m) - d/dy(t,/m)], 
where m is the map factor (see Haltiner and Williams, 1980): 
m(~) = (1 + sin 60°)/(1 + sin), 


and t, and ty are in the polar stereographic reference frame. Wind stress time series 


from four (in a square configuration) grid points or three (in a triangle) may be used to 


ie 


compute finite difference approximations to the derivatives in the equation for V x t 
above. The geographic center of the current meter array used in this study is positioned 
with respect to the FNOC grid such that use of the square scheme to compute wind 
stress curl requires that a grid point over land be used in the computation. Use of the 
triangular scheme allows exclusion of the grid point over land from the computation. 
Wind stress curl time series using both of these schemes were analyzed to see if 
the larger truncation error associated with the triangular scheme was significant. The 
time series derived from the square (F-G-J-K, Figure 2.1) and triangular (F-J-K, Figure 
2.1) schemes look virtually identical, Figure A.1. The autospectra of the V x t time 
series derived from the square and triangular schemes are very similar with only minor 
differences detectable, Figure A.2. The time series of V x t from these schemes are 
coherent for all periods less than 90 days, Figure A.3. Hence, the effect of additional 


truncation error of the triangular scheme 1s not readily apparent. 
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Figure A.1 Time series of wind stress curl derived from 
Square (top) and triangular (bottom) schemes. 
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Figure A.2. Autospectra of wind stress curl computed with 
Square (a) and triangular (b) schemes. 
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Figure A.3. Coherence and phase between square and 
triangular schemes. 
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